Consistent testing for a constant copula under strong mixing 
based on the tapered block multiplier technique 
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^ ■ Abstract 



Considering multivariate strongly mixing time series, nonparametric tests for a constant 
copula with specified or unspecified change point (candidate) are derived; the tests are con- 
sistent against general alternatives. A tapered block multiplier technique based on serially 
dependent multiplier random variables is provided to estimate p- values of the test statistics. 
Size and power of the tests in finite samples are evaluated with Monte Carlo simulations. 
The block multiplier technique might have several other applications for statistical inference 
on copulas of serially dependent data. 
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1. Introduction 



Over the last decade, copulas have become a standard tool in modern risk management. 
The copula of a continuous random vector is a function which uniquely determines the 
dependence structure linking the marginal distribution functions. Copulas play a pivotal 
role for, e.g., measuring multivariate association [see 146], pricing multivariate options [see 
49l | and allocating financial assets [see . The latter two references emphasize that time 
variation of copulas possesses an important impact on financial engineering applications. 

Evidence for time-varying dependence structures can indirectly be drawn from functionals 
of the copula, e.g., Spearman's p, as suggested by Gaifeer et al. [20|_and Wied et al. 52|. 
Investigating time variation of the copula itself, Busetti and Harvey [9fl consider a nonpara- 
metric quantile-based test for a pointwise constant copula. Semiparametric tests for time 
variation of the parameter within a prespecified parametric copula model are proposed by 
Dias and Embrechts fl5| and Giacomini et al. (23 |. Guegan and Zhang 24 1 combine tests 
for constancy of the copula (on a given set of vectors on its domain), the copula family, 
and the parameter. All the latter references on time variation of the copula are based on 
the assumption of i.i.d. observations. With respect to financial time-series, this assumption 
may be approximated by the estimation of a GARCH model and by using the residuals 
obtained after GARCH filtration. However, the effect of replacing unobserved innovations 
by estimated residuals has to be taken into account. Therefore, specific techniques for 
residuals are required (cf. 01) and exploring this approach, Remillard investigates a 
nonparametric change point test for the copula of residuals in stochastic volatility models. 

In the present paper we go beyond these approaches and work in a purely nonparametric 
setting and under the allowance of general serial dependence of a multivariate time series 
measured by their al pha -mixing coefficients. In similar settings, Fermanian and Scaillet 
(l9| . Doukhan et al. [17] and Biicher and Volgushev Q consider the nonparametric esti- 
mation of copulas. These references form the basis for our new tests for a constant copula 
under strong mixing mixing assumptions. We introduce nonparametric Cramer- von Mises-, 
Kuiper-, and Kolmogorov-Smirnov tests which assess constancy of the copula on its entire 
domain. In consequence, they are consistent under general alternatives. Depending on the 
object of investigation, tests with a specified or unspecified change point (candidate) are 
introduced. Whereas the former setting requires a hypothesis on the change point location, 
it allows us to relax the assumption of strictly stationary univariate processes. P- values of 
the tests are estimated based on a generalization of the multiplier bootstrap technique in- 
troduced in Remillard and Scaillet [40| to the case of strongly mixing time series. This new 
technique can be used to generalize several statistical inference methods from the i.i.d. case 
to the serial dependent case. The idea of our approach is comparable to block bootstrap 
methods with the difference that, instead of sampling blocks with replacement, we gener- 
ate blocks of serially dependent multiplier random variables. For a general introduction to 
the latter idea, we refer to Buhlmann (?J and Paparoditis and Politis (3l| . Independently 
from our work, a current working paper of van Kampen and Wied |5l| investigates similar 
Kolmogorov-Smirnov tests for constancy of the copula with unspecified change point based 
on a bootstrap technique introduced in Inoue [26]. 

This paper is organized as follows. In Section [21 we discuss weak convergence of the em- 
pirical copula process under strong mixing. We derive a tapered block multiplier bootstrap 
technique for inference on the weak limit and assess it in finite samples. Tests for a con- 
stant copula with specified or unspecified change point (candidate) which are relying on 
this technique are established in Section O Section 0] concludes the paper. All the proofs 
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are deferred to 



Appendix A 



2. Nonparametric inference based on serially dependent observations 

As a basis for the tests introduced in Section[3j we briefly recapitulate the results of Doukhan 
et al. (l~7| . Segers 47 1 and Biicher and Volgushev on the asymptotic behavior of the em- 



pirical copula process under nonrestrictive smoothness assumptions in the case of strongly 
mixing observations. The main result of this section is a multiplier-based resampling method 
for this particular setting. We establish its asymptotic behavior and investigate the perfor- 
mance in finite samples. 

2.1. Asymptotic theory 

Consider a vector- valued process (Xj)j 6 z with Xj = (Xj t i, . . . , -Xj,<f) taking values in W 1 . 
Let Fi be the distribution function of Xji for all j £ Z, i = 1,. . . ,d and let F be the 
joint distribution of Xj for all j £ Z. Assume that all marginal distribution functions are 
continuous. Then, according to Sklar's Theorem [48], there exists a unique copula C such 
that F(x\, . . . , Xd) = C(F\(xi), . . . , Fd(xd)) for all (xi, . . . , Xd) £ K d . The a-fields generated 
by Xj, j < t, and Xj, j > t, are denoted by Ft = <r{Xj, j < t} and F 1 = <j{Xj,j > t}, 
respectively. We define 

a(F s ,F s+r ) = sup \P(AnB)-P(A)P(B)\. 

The strong- (or a-) mixing coefficient ax corresponding to the process (Xj)j g z is given by 
ax(r) = sup s>0 a(F s ,F s+r ). The process (Xj-)j-gz is sa id to be strongly mixing if ax(r) — > 
for r — > oo. This type of weak dependence covers a broad range of time-series models. 



Consider the following examples, cf. Doukhan [16| and Carrasco and Chen 

Example 1. i) AR(1) processes (Xj)j £ i given by 

Xj = /3X,-_i + Ej, 

where (£j)j£% is a sequence of independent and identically distributed continuous innovations 
with mean zero. For |/3| < 1, the process is strictly stationary and strongly mixing with 
exponential decay of «x(r). 

ii) GARCH(1, 1) processes (Xj)j £ z, 

Xj = o-jEj, = oj + /3o"j_i + aef-i, (1) 

where {sj)j^% is a sequence of independent and identically distributed continuous innova- 
tions, independent of (Tq, with mean zero and variance one. For a + (3 < 1, the process is 
strictly stationary and strongly mixing with exponential decay o/ax(r). 

Hi) For multivariate analogues of i) and ii) see Section \2. 3\ below. 

Let Xi, . . . ,X n denote a sample from (Xj)j e z- A simple nonparametric estimator for the 
copula C is given by the empirical copula which is first considered by Riischendorf |44j and 



Deheuvels [14J. Depending on whether the marginal distribution functions are assumed to 



be known or unknown, we define, for u £ [0, l] d , 



^ n 1 n 

C »M ■■= -EV,<u}, C n (u) := - £l (2) 



n ^ l u i< u }' 
i=i j=i 
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where Uj = (Ujx, . . . , Uj^) and Uj = (Uj i, . . . ,Uj 4) with observations U»i = Fi(Xj^) 

and pseudo-observations Uj^ = Fi(Xjj) for all j = 1, . . . ,n and i = l,...,d, and where 

Fi{x) = n -1 Ysj=i ^-{Xj i<x} f° r au x ^ ^- Unless otherwise noted, the marginal distribution 

functions are assumed to be unknown and the empirical copula C n is used. In addition 
to the practical relevance of this assumption, Genest and Segers [22] prove that pseudo- 
observations Uj permit more efficient inference on the copula than observations Uj for a 
broad class of copulas. 

Doukhan et al. fl~7| investigate dependent observations and establish the asymptotic behav- 
ior of the empirical copula process, defined by \/n{C n — C}, assuming the copula to possess 
continuous partial derivatives on [0, l] d . Segers 47J] points out that many popular families 
of copulas (e.g., the Gaussian, Clayton, and Gumbel-Hougaard families) do not satisfy the 
assumption of continuous first partial derivatives on [0, l] d . He establishes the asymptotic 
behavior of the empirical copula for serially independent observations under the weaker 
condition 

DjC(u) exists and is continuous on |u G [0, l] d \ui G (0, 1) \ for all i = 1, . . . , d. (3) 
Under Condition ([3]), the partial derivatives' domain can be extended to u G [0, l] d by 

lim^o C(u+fee ^~ C(u) for all u G [0, l] d , < u t < 1, 
limsup^o C(u + fce ' } for all u G [0, l] d , m = 0, (4) 

k limsup^Q C(u) ~ C /t (u ~ fce ' ) for all u G [0, m = 1, 

and for all i = 1, . . . , d, where denotes the ith column of a d x d identity matrix. The 
following generalization of the results in 17] and 47| is a consequence of Theorem 2.4 in 
Biicher and Volgushev jfj, see Corollary 2.5 in that reference. 



D l C(u) 



Theorem 1. Consider observations Xi, . . . ,X n , drawn from a strictly stationary process 
(Xj)j g z satisfying the strong mixing condition ax(r) = 0(r~ a ) for some a > 1. Then 

B C ,n = Vn~(C n ~ C) B c 

in the metric space space of uniformly bounded functions on [0, l] d equipped with the uniform 
metric (£°°([0,l] d ), \\ ■ ||oo)- Here, Be* denotes a centered tight Gaussian process on [0, l] d 
with covariance function 

7 (u,v) = Cov(B c (u),B c (v)) =^Cov(l {Uo < u} ,l {u .< v} ) for all u, v G [0, l] d . (5) 

Moreover, if C satisfies Condition |3|) ; then 

Gc,n = Vn(C n - C) G c 

in £°°([0, l] d ), where Gc represents a Gaussian process given by 

d 

G c (u) = B c (u) - AC(u)B c (u«) for all u G [0, l] d . (6) 

Here, denotes the vector where all coordinates, except the ith coordinate of u, are 

replaced by 1. 

Notice that the covariance structure as given in Equation iJS]) depends on the entire process 
(Xj)j g z in case it is not serially uncorrelated. 
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2.2. Resampling techniques 



In this Section, we introduce two bootstrap techniques for the empirical copula process 
which are applicable in the case of strongly mixing observations. We begin with the (mov- 
ing) block bootstrap, which serves as a benchmark in the finite sample assessment. Sub- 
sequently, we derive the main result of this Section about the asymptotic consistency of a 
generalized multiplier bootstrap technique. 



2.2.1. The block bootstrap 

Fermanian et al. [18] investigate the empirical copula process for independent and identi- 
cally distributed observations Xi, . . . , X n and prove consistency of the nonparametric boot- 
strap method which is based on sampling with replacement from Xi, . . . , X n . We denote a 
bootstrap sample by X;f , . . . , X^ and define 

d n («) : = - E !{0f <»} fOT a11 U G t°> Vf* H ■= ~ E \XS 

j=l k=l 

for all j = 1, . . . , n and i = 1, . . . , d. Notice that the bootstrap empirical copula can equiva- 
lently be expressed based on multinomially (n,n _1 , . . . ,n~~ l ) distributed random variables 
W = (W 1 ,...,W n ): 

1 n n 

d n(") : = -E^ 1 {U f <u} fOT a11 U 6 M*. U% ■■=-Y, W ^{X k , i <X j , i } 
3=1 k=l 

for all j = 1, ... ,n and i = 1, ... ,d. It is shown in (l8l . Q that the corresponding boot- 
strap empirical copula process converges weakly conditional on Xi, . . . ,X n in probability in 
(£°°([0, l] d ), || • \\oo), notationally 

= yfcffi - C n ) A Go. (7) 

Here, weak convergence conditional on the data in probability is understood in the Hoffmann- 
J0rgensen sense as c 
two conditions hold: 



J0rgensen sense as defined in Kosorok [29(|, i.e., G^ n — > Gc if and only if the following 



sup 

heBLi(£°°([0,l] d )) 



E w h(G^ n ) - Eh (G c ) -^O, (8) 



E w h{G^ n )* - E w h{G^ n \ -A for all h G BLi(*°°([0, l} d )). (9) 

p 

Here, — > denotes convergence in outer probability and Ey/ denotes expectation with 
respect to W conditional on Xi, . . . ,X n . Furthermore, BL±(£°°([0, l] d )) denotes the set 

/ . f»([ , 1]d) ^ R i < 1; \ f{f3) _ /(7) | < | !/3 _ 7 || Qo for all 7) ^ G ^oc ([0j 

of all Lipschitz-continuous functions bounded by 1 with Lipschitz-constant not exceeding 
1 and the asterisks in ([9]) denote measurable majorants and minorants with respect to the 
joint data (i.e., Xi, . . . , X n and W). Weak convergence conditional on Xi, . . . , X n almost 
surely is defined analogously by replacing outer probability convergence in ([8]) and ([9]) by 
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outer almost sure convergence. Due to the lack of a general continuous mapping Theorem 
and an easy functional delta method, see 50, 2^|, we do not consider the outer almost sure 



version in this paper. 

Whereas the bootstrap is consistent for i.i.d. samples, consistency generally fails for serially 
dependent samples. Therefore, a block bootstrap method is proposed by Kiinsch (3oT |. Given 
the sample Xi, . . . ,X n , the block bootstrap method requires blocks of size Is = lB( n ), 
lB( n ) — )• oo as n — )• oo and Ib(ji) = o(n), consisting of consecutive observations 

B h,l B = {Xh+i, . . . ,X h +i B }, for all h = 0, . . . ,n- Ib- 

We assume n = Mb (otherwise the last block is truncated) and simulate H = {Hi, . . . , H^) 
independent and uniformly distributed random variables on {0, ...,n — Ib}- The block 
bootstrap sample is given by the observations of the k blocks Bh x i b , ■ ■ ■ ,Bjj i B , i.e., 

Xffi+ij . . . ,X.H 1+iB ,X.H- 2+ i, . . . ,X. H2+ i g , ,X# fe +i, • • • ,X Hk+ i B . 

Denote the block bootstrap empirical copula based on this sample by C^(u). Its asymptotic 
behavior can be established by means of an asymptotic result on the block bootstrap for 
general <i-dimensional distribution functions established by Buhlmann 0, Theorem 3.1], see 
also Example 2.10 in Biicher and Volgushev 0|. 

Theorem 2. Consider observations Xi, . . . ,X n , drawn from a strictly stationary process 
(Xj)jez satisfying YlT=i( r + l) 16(d+1) v / ^xH < 00 • Assume that l B (n) = 0{n l / 2 ~ £ ) for 
< e < 1/2. If C satisfies Condition then the block bootstrap empirical copula process 
converges weakly conditional on Xi, . . . , X n in probability in (£°°([0, l] d ), || ■ ||oo) : 

Gg n (u) = V^{d*(u) - d n (u)} A Go(u). 

A brief remark on the condition on the mixing rate is in order here. The summability 
conditions forces ax(r) to be of order 0(r~ 96_£ ) for dimension d = 2, which is far away 
from the (sharp) rate 0(r -1 ~ e ) needed in Theorem [TJ This discrepancy is due to the 
fact that the literature does not provide stronger results on the consistency of the block 
bootstrap for the d-variate empirical process except the ones in [7|, at least to the best of 
our knowledge. Exploiting more recent techniques (which are beyond of the scope of this 
paper) we believe that it is possible to get better rates which are comparable to those of the 



non-bootstrap version, see also [35j| and 37] . The proof of Theorem [2] being based on the 
functional delta method would easily transfer these rates to the rank-based copula setting 
considered in the present paper. Also note that for applications the message is not too 
bad: most time series models have exponentially decreasing alpha mixing coefficients, see 
Example [T] or more precisely the examples in 0| . 

2.2.2. The multiplier bootstrap 

A process related to the bootstrap empirical copula process defined in ([7]) can be formulated 
if both the assumption of multinomially distributed random variables is dropped and the 
marginal distribution functions are left unaltered during the resampling procedure. The 
resulting bootstrap scheme is known for the i.i.d. context as the multiplier bootstrap or 



the multiplier method, see [45l . |4Q, |4|, |47| . In the present section we briefly summarize this 



concept and extend it to the serially dependent case. 
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Consider i.i.d. multiplier random variables £1, . . . , £ n with mean and variance 1, additionally 
satisfying ||£j||2,i := V^G^I > x)dx < oo for all j = 1, ... , n (where the last condition 
is slightly stronger than that of a finite second moment). Replacing the multinomial mul- 
tiplier random variables Wi, . . . , W n by . . . , where £ = n" 1 Y^j=i £ji (ensuring 
realizations having arithmetic mean one) yields the multiplier (empirical copula) process 
which converges weakly conditional on Xi, . . . ,X n in probability in (£°°([0, l] d V, || • ||oo) in 
the i.i.d. situation, see Biicher and Dette [4j and more precisely in Biicher et al. p, Theorem 
2.3]: 



n < 



„ £ {U,<u} 



n * — ' £ 



C„(u) 



u 



For general considerations of multiplier empirical processes, we refer to the monographs of 
van der Vaart and Wellner 50j| and Kosorok [29j. The process is introduced by Scaillet 
45| in a bivariate context, a general multivariate version and its unconditional weak con- 
vergence are investigated by Remillard and Scaillet [40] and Segers [47]. In order to get 



approximations of the limiting process Gc these authors propose to estimate the partial 
derivatives of the copula in ([6]) by some estimator DiC n (see, e.g., Example [2] below) and 
define 



"C,n 



U 



which converges to Gc conditionally on the data in probability, see (3, 0|. Here, the estima- 
tor DiC n is supposed to satisfy the following two assumptions, which will also be necessary 
in the serially dependent situation: 

CI There exists a constant K such that ||-DjC n ||oo < K for all n G N. 
C2 For all 5 G (0, 1/2) one has 



sup 

u6[0,l] d :u 4 e [<5,l-<5] 



AC„(u) - AC(u; 



0. 



Example 2. is easily seen that finite differencing yields a simple nonparametric estima- 
tor for the first order partial derivatives DjC(u) which satisfies CI and C2. More precisely, 
we define 



AC(u) 



C„(u+^)-C„(u-/,e0 for aU u £ [Qj < n . < ! _ 

C„(u+2/» ei ) for aU u G ^ ^rf q < u . < ^ 

Cn(u)-C n (u-2feei) f _ ^ n< ^ fn lid 
2ft 



for all u G [0, 1 - /i< ^ < 1, 



(10) 



where h = h n such that inf n h n y/n > and where denotes the ith column of the 
d x d identity matrix [see 47. 3|. 

Biicher and Dette (3| find that the multiplier technique yields more precise results than 
the nonparametric bootstrap in mean as well as in mean squared error when estimating 
the asymptotic covariance of the empirical copula process in the i.i.d. context. Motivated 
by the fact that this technique is inconsistent when applied to serially dependent samples, 



G 



we derive a generalization of the multiplier technique in the following. Inoue 26] develops 
a block multiplier process for general distribution functions based on dependent data in 
which the same multiplier random variable is used for a block of observations and the 
composition of blocks remains unaltered throughout the procedure. We consider a technique 
in which the composition of blocks is not fixed. More precisely, a tapered block multiplier 
(empirical copula) process is introduced based on the work of Biihlmann [7, Chapter 3.3] 



and Paparoditis and Politis 3lJ: The main idea is to consider a sample £i n , . . . , ^ n ^ n from a 
process (£j,n)jez of serially dependent tapered block multiplier random variables, satisfying: 

Al (£j,n)jez is independent of the observation process (Xj)j S z- 

A2 (£j n)jeZ i s a positive c • I (n) -dependent process, i.e., for fixed j £ Z, £j n is independent 
of ij+h,n f° r all \h\ > c • l(n), where c is a constant and l{n) — > oo as n — > oo while 
l{n) = o{n). 

A3 (Cj,n)jeZ is strictly stationary. For all j, h E Z, assume n] = A 4 > 0, Cov[^.„, £j+h,n] = 
H 2 v(h/l(n)) and f is a bounded function symmetric about zero; without loss of generality, 
we consider [i = 1 and v(0) = 1. All central moments of £j n are supposed to be bounded. 

The following theorem is the main result of this section. Weak convergence of the tapered 
block multiplier process conditional on a sample Xi, . . . , X n is established. Regarding the 
strong assumptions on the mixing rate the remark after Theorem [2] holds true here as well. 

Theorem 3. Consider observations Xi,...,X n drawn from a strictly stationary process 
(Xj)j g z satisfying X^i( r + l) c \/ a x(^) < oo, where c = max{8d + 12, [2/e\ + 1}. Let the 
tapered block multiplier process (£j, n )je% satisfy Al, A2, A3 with block length l(n) — > oo, 
where l(n) = 0{n 1/2 ~ e ) for < e < 1/2. Then, 



^(u)=^fif:^i { o i <„ } -^(u) 



u) in l°°([0,l] d ), 



where £ ra = n z^?=i£i,n- Moreover, if ([3]) holds and if DiC n satisfies conditions CI and 
C2, then 



n (u) = Bg n (u) -^A^(u)B^ (u«) A G c (u) ? n £°°([0, 1] 



Remark 1. The multiplier random variables can as well be assumed to be centered around 
zero [cf. 29, Proof of Theorem 2.6]. Define £j n := £j >n — /i. Then 



l C. 



•'W = Tn g ( «» " ® 'PW " ^ g (If " ^ " |S "" (U) 

/or a// u G [0, l] d . 77us is an asymptotically equivalent form of the above tapered block 
multiplier process: 



sup 

[o,i] d 



sup 

[0,1] d 



(C„-l) B#„(u) 







almost surely, since B^ n (u) tends to a tight centered Gaussian limit, unconditionally. The 
assumption of centered multiplier random variables is abbreviated as A3b in the following. 
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In practice, given observations Xi, . . . , X n of a strictly stationary process (Xj)j 6 z satisfying 
the assumptions of Theorem [3l approximating Gc* requires the following three steps: 

1. Estimate the partial derivatives of C by an estimator satisfying conditions CI and 
C2, for instance by the estimator given in Example [2] above. 

is) (s) 

2. For s = 1, . . . , S with S G N, simulate samples Ci nj • • • ; £n,n from a tapered block 
multiplier process (£j in )jez satisfying Al, A2, A3. For each s G {1, . . . , S 1 }, calculate 

1 n /f {s) \ 

s 5 s) h = -i E |S " 1 Mi<u } for a11 u G [°. ^ (") 

3. For each s G {1, . . . , S}, calculate 

egW(u) = B#W(u) - (u«) for all u G [0, (12) 

i=l 

The sample {G c f^} s= i j ....5' is an approximate sample of Gc- 

There are numerous ways to define tapered block multiplier processes (£j,n)jez satisfying 
the assumptions Al, A2, A3 made in Theorem [3l In the remaining part of this section, 
a basic version having uniform weights and a refined version with triangular weights are 
investigated and compared. 

Example 3. A simple form of the tapered block multiplier random variables can be defined 
based on moving average processes. Consider the function K\ which assigns uniform weights 
given by 

Kl (h) : =fw=I f° r aU \ h \ < Kn) 
I else. 

Note that K\ is a discrete kernel, i.e., it is symmetric about zero and X^heZ ^iC 1 ) = 1- ^ e 
tapered block multiplier process is defined by 

oo 

= Cj,n = E K ii h ) w j+h for all j G Z, (13) 

h=— oo 

where (wj)j^z * s an independent and identically distributed sequence of, e.g., Gamma(q,q) 
random variables with q := l/[2Z(n) — 1]. The expectation of £j is then given by E[£j] = 1, its 
variance by Var[^j] = 1 for all j G Z. For all j G Z and \h\ < 2l(n) — 1, direct calculations 
further yield the covariance function Cav(£j,£ J \f/ l ) = {21 (n) — 1 — |/i|}/{2Z(n) — 1} which 
linearly decreases as h increases in absolute value. The resulting sequence (Cj)jeZ satisfies 
Al, A2, and A3. Exploring Remark{l\ tapered block multiplier random variables can as 
well be defined based on sequences (wj)j£z °fi e -9-, Rademacher-type random variables Wj 
characterized by P(u>j = —1/y/q) = P(wj = l/^/q) = 0.5 or Normal random variables 
Wj ~ N(0, 1/y/q). In either one of these two cases, the resulting sequence (Cj)jeZ satisfies 
Al, A2, and A3b. Figure [7] shows the kernel function K\ and simulated trajectories of 
Rademacher-type tapered block multiplier random variables. 
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Figure 1: Tapered block multiplier Monte Carlo simulation. Kernel function K\{h) (left) and 
simulated trajectories of Rademacher-type tapered block multiplier random variables £1,... , £100 
(right) with block length l(n) = 3 (solid line) and l(n) = 6 (dashed line), respectively. 
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Figure 2: Tapered block multiplier Monte Carlo simulation. Kernel function K2Q1) (left) and 
simulated trajectories of Rademacher-type tapered block multiplier random variables £1,... ,£100 
(right) with block length l(n) = 3 (solid line) and l(n) = 6 (dashed line), respectively. 



Example 4. Following Biihlmann f^J, let us define the kernel function by 

K 2 (h) := max{0,{l - \h\/l(n)}/l(n)} 

for allh 6 Z. The tapered multiplier process (Cj)jeZ follows Equation t!3\) . where (vjj)j^z * s 
an independent and identically distributed sequence of Gamma(q,q) random variables with 
q = 2/{3l(n)} + l/{3^(n) 3 }. The expectation of £j is given by 



l(n) . 



l(n) 



1. 



For the variance, direct calculations yield 

l(n) 



VarU 



l(n 



{l(n)-hY 



h=l 



l(n) 



Var[w, 



+ 



3l(n) 3l(n) 3 



Var\w 



for all j G Z. For any j £ Z and \h\ < 2l(n) — 1, the covariance function Cov(^j, £j+h) 
can be described by a parabola centered at zero and opening downward [for details, see l& 
Section 6.2]. The resulting sequence (Cj)jeZ satisfies Al, A2, and A3. Figure^ provides an 
illustration of the kernel function K2 as well as simulated trajectories of Rademacher-type 
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tapered block multiplier random variables. In this setting, (£j)jez satisfies Al, A2, and 
A3b. Notice the smoothing which is driven by the choice of kernel function and the block 
length l(n). This effect can be further explored using more sophisticated kernel functions, 
e.g., with bell-shape; this is left for further research. 



2.3. Finite sample behavior 

In the present section we investigate and compare the finite sample properties of the (mov- 
ing) block bootstrap and the tapered block multiplier technique by means of a simula- 
tion study. For that purpose we use both resampling techniques in order to estimate the 
(co)variances of the empirical copula process. This study complements the one of Bucher 
and Dette [4] and Bucher 0| on bootstrap approximations for the empirical copula process 
in the i.i.d setting. 

To be more precise, in Tables [U [2] and [3] we demonstrate simulation results (for n = 100 and 
n = 200) on the estimation of the theoretical covariance Cov(Gc(u), Gc(v)) calculated at 
the points u = v G {(1/3, 1/3), (2/3, 1/3), (1/3, 2/3), (2/3, 2/3)} in basically six different 
models: serial dependence features arise from (bivariate) i.i.d., AR(1) or GARCH(1, 1) time 
series models, while the copulas linking the marginals of the innovations are taken from the 
Gumbel-Hougaard or the Clayton family. The technical details of the execution of the 
study are given below, we start with the discussion of the results. 

The results based on independent and identically distributed samples indicate that the ta- 
pered block multiplier outperforms the block bootstrap in mean and MSE of estimation. 
However, applying the resampling methods of the present paper in an i.i.d. setting comes 
at the price of an slight increased mean squared error in comparison to the multiplier or 
(block) bootstrap with block length I = 1 as investigated in [4|. Hence, we suggest to test 
serial independence of continuous multivariate time-series as introduced by Kojadinovic and 
Yan [28] to investigate which method is appropriate. In the case of serially dependent ob- 
servations, the results indicate that the tapered block multiplier yields more precise results 
in mean and mean squared error than the block bootstrap (which tends to overestimate) 
for the considered choices of the temporal dependence structure, the kernel function and 
the copula. This finding coincides with the results in ji| . Regarding the choice of the kernel 
function, mean results for K\ and k 2 are similar, whereas k 2 yields slig htly better results in 
mean squared error. Additional MC simulations are given in Ruppert [43]: if the multiplier 
or bootstrap methods for independent observations are incorrectly applied to dependent 
observations, i.e., Ib = Im = 1) then their results do not reflect the changed structure 
adequately. Results based on Normal, Gamma, and Rademacher-type sequences (vjj)j£z 
indicate that different distributions used to simulate the multiplier random variables lead 
to similar results. To ease comparison of the next section with the work of Remillard and 



Scaillet [40j], we decided to state the results for normal multiplier random variables here. 

Finally, for the sake of completeness we state the technical details underlying the simulation. 
First of all, the details regarding the models are as follows. 

• The bivariate Clayton and Gumbel-Hougaard copulas are given by 
C$ l {u l ,u 2 ) = (V e + u 2 e , 6>0, 



Cf>i,u 2 ) =exp( - {-ln( Ul )Y + {-ln(u 2 )y 
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Table 1: Mean and MSE (xlO 4 ) Monte Carlo results. I.i.d. and AR(1) settings, sample 
size n = 100 and 1,000 Monte Carlo replications. For each replication, we perform S = 2,000 
tapered block multiplier (Mi) repetitions with Normal multiplier random variables, kernel function 
Ki, i = 1,2, block length Im = 3, and block bootstrap (B) repetitions with block length Ib = 5. 



(ui,u 2 ) 




(1/3,1/3) 




(1/3,2/3) 




(2/3,1/3) 




(2/3,2/3) 








Mean 




MSE 


Mean 




MSE 


Mean 




MSE 


Mean 




MSE 


i.i.d. setting 


























Clayton 


True 


0.0486 






0.0338 






0.0338 






0.0508 






(9 = 1) 


Approx. 


0.0487 






0.0338 






0.0338 






0.0508 








M 2 


0.0496 


1 


.6323 


0.0344 


1 


.2449 


0.0345 


1 


.3456 


0.0528 


1 


.6220 




Mi 


o n/iG/i 


1 
1 


.7949 




1 
1 


.2934 


G HQ/1 Q 


1 
1 


.4128 


U.U0Z4 


1 
1 


.7910 




B 


0.0599 


2 


.8286 


0.0432 


2 


.4103 


0.0429 


2 


.1375 


0.0643 


3 


.1538 


Clayton 


True 


0.0254 






0.0042 






0.0042 






0.0389 






(9 = 4) 


Approx. 


0.0255 






0.0042 






0.0042 






0.0390 








M 2 


0.0259 





.9785 


0.0051 





.3715 


0.0048 





.3662 


0.0407 


1 


.6324 




M 1 


09^7 


1 
1 


.0104 


u.uuou 


n 
u 


.3656 


U.UU4o 


n 
u 


.3673 


A 0/10/1 
U.U4U4 


1 
1 


.7048 




B 


0.0383 


2 


,5207 


0.0100 





.7222 


0.0097 





.6662 


0.0533 


3 


.7255 


Gumbel 


True 


0.0493 






0.0336 






0.0336 






0.0484 






(9 = 1.5) 


Approx. 


0.0493 






0.0335 






0.0335 






0.0485 








M 2 


0.0514 


1 


.3914 


0.0346 


1 


.2657 


0.0340 


1 


.2583 


0.0497 


1 


.4946 




Mi 


U.UOlU 


1 
1 


.5530 


H GQ/l A 
U.Uo44 


1 
1 


.3390 


U.Uooo 


1 
1 


.3275 


n O/i 


1 
1 


.6948 




B 


0.0616 


2 


.8334 


0.0429 


2 


.1366 


0.0432 


2 


.2273 


0.0620 


3 


.2134 


Gumbel 


True 


0.0336 






0.0058 






0.0058 






0.0293 






(9 = 3) 


Approx. 


0.0335 






0.0058 






0.0058 






0.0294 








M 2 


o hq^p; 

U.UoOO 


1 
1 


.1359 


U.UU04 


u 


.4427 


U.UUDo 


u 


.3851 


n 0^07 


u. 


.9819 




Mi 


0.0353 


1 


.1991 


0.0063 





.4409 


0.0062 





.3832 


0.0306 


i 


.0261 




B 


0.0470 


2 


. ( ooo 


0.0120 





.ooyu 


0.0122 







0.0437 


2 


• if ooo 


AR(1) setting with 


P = 0.5 
























Clayton 


Approx. 


u.uoyy 






U.U4Uo 






n G/1 no 
u.U4uy 






u.uozy 






(9 = 1) 


M 2 


0.0602 


3 


.1797 


0.0394 


2 


.4919 


0.0398 


2 


.5903 


0.0625 


2 


.6297 




Mi 


0.0598 


3 


■ OOVO 


0.0391 


2 


.uuyu 


0.0396 


2 


7zL1 n 


0.0620 


2. 


.yozu 




B 


0.0699 


3 


.4783 


0.0496 


2 


.9893 


0.0492 


2 


.8473 


0.0761 


4 


.0660 


Clayton 


Approx. 


0.0329 






0.0064 






0.0064 






0.0432 






(9 = 4) 


M 2 


0.0347 


2 


.0017 


0.0071 





.6040 


0.0072 





.6666 


0.0460 


2. 


.8656 




Mi 


0.0344 


2 


.0672 


0.0071 





.5869 


0.0072 





.6583 


0.0458 


3 


.0571 




B 


0.0472 


3 


.5289 


0.0132 


1 


.0990 


0.0128 





.9658 


0.0585 


4 


.2965 


Gumbel 


Approx. 


0.0617 






0.0408 






0.0409 






0.0605 






(9 = 1.5) 


M 2 


0.0631 


3 


.0587 


0.0406 


2 


.7512 


0.0398 


2 


.4187 


0.0600 


3 


.0433 




Mi 


0.0626 


3 


.4252 


0.0402 


2 


.8739 


0.0395 


2 


.4920 


0.0594 


3 


.3824 




B 


0.0735 


3 


.7410 


0.0501 


3 


.2025 


0.0502 


3 


.1150 


0.0735 


4 


.0521 


Gumbel 


Approx. 


0.0385 






0.0061 






0.0061 






0.0345 






(9 = 3) 


M 2 


0.0425 


2 


.5441 


0.0069 





.5841 


0.0070 





.5796 


0.0375 


2 


.1518 




Mi 


0.0422 


2 


.7239 


0.0069 





.5706 


0.0070 





.5836 


0.0373 


2. 


.2381 




B 


0.0535 


3 


.8803 


0.0127 


1 


.0720 


0.0125 


1 


.0526 


0.0507 


4 


.1894 
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Table 2: Mean and MSE (xlO 4 ) Monte Carlo results. I.i.d. and AR(1) settings, sample 
size n = 200 and 1,000 Monte Carlo replications. For each replication, we perform S = 2,000 
tapered block multiplier (Mj) repetitions with Normal multiplier random variables, kernel function 
Ki, i = 1,2, block length Im = 4, and block bootstrap (B) repetitions with block length Ib = 7. 



(ui,u 2 ) 




(1/3,1/3) 




(1/3,2/3) 




(2/3,1/3) 




(2/3,2/3) 








Mean 




MSE 


Mean 




MSE 


Mean 




MSE 


Mean 




MSE 


i.i.d. setting 


























Clayton 


True 


0.0486 






0.0338 






0.0338 






0.0508 






(9 = 1) 


Approx. 


0.0487 






0.0338 






0.0338 






0.0508 








M 2 


0.0494 


1 


.2579 


0.0346 





.8432 


0.0343 





.8423 


0.0522 


1 


.0987 




Mi 


o n/ion 
U.U4yu 


1 
1 


.3749 


U.Uo40 


n 
u 


.8880 


U.Uo41 


n 
u 


.8719 


O or;i g 


1 
1 


.2074 




B 


0.0562 


1 


.7155 


0.0402 


i 


.2154 


0.0395 


i 


.2279 


0.0593 


1 


.8137 


Clayton 


True 


0.0254 






0.0042 






0.0042 






0.0389 






(9 = 4) 


Approx. 


0.0255 






0.0042 






0.0042 






0.0390 








M 2 


0.0261 





.5811 


0.0044 





.1889 


0.0048 





.2027 


0.0390 





.9932 




Mi 


o nof^n 
u.uzou 


n 
u 


.6024 


U.UU44 


n 
u 


.1876 


U.UU4o 


u 


.2008 


0^88 


1 
1 


.0702 




B 


0.0347 


i 


.6381 


0.0076 





.3197 


0.0073 





.2908 


0.0489 


2 


.1425 


Gumbel 


True 


0.0493 






0.0336 






0.0336 






0.0484 






(9 = 1.5) 


Approx. 


0.0493 






0.0335 






0.0335 






0.0485 








M 2 


0.0512 


i 


.1513 


0.0347 





.8157 


0.0347 





.8113 


0.0504 


1 


.1697 




Mi 


o gp;i 1 


i 
i 


.2909 




u 


.8653 


U.Uo4D 


n 
u 


.8692 


O O^OQ 


1 

1 


.2811 




B 


0.0577 


i 


.7178 


0.0402 


i 


.2233 


0.0403 


i 


.1701 


0.0574 


1. 


.8245 


Gumbel 


True 


0.0336 






0.0058 






0.0058 






0.0293 






(9 = 3) 


Approx. 


0.0335 






0.0058 






0.0058 






0.0294 








M 2 


o oqai 


n 
u 


.8340 


o c\c\c\i 

U.UUO i 


u 


.2577 


n oo^Q 


u 


.2505 


O OQOO 


u. 


.8678 




Mi 


0.0359 





.9107 


0.0067 





.2576 


0.0062 





.2444 


0.0318 





.9078 






0.0435 


1 




0.0095 





7 


0.0094 





.oou^t 


0.0388 


1 




AR(1) setting with 


P = 0.5 
























Clayton 


Approx. 


u.uoyy 






U.U4Uo 






n O/i on 
u.U4uy 






O OfiOG 






(9 = 1) 


M 2 


0.0615 


2 


.3213 


0.0413 


1 


.8999 


0.0414 


1 


.9235 


0.0646 


2 


.2534 




Mi 


0.0608 


2 


c. c; c. o 


0.0408 


1 




0.0410 


2 




0.0638 


2 








0.0676 


2 


.6206 


0.0468 


1 


.9802 


0.0472 


2 


.0278 


0.0715 


2 


.5494 


Clayton 


Approx. 


0.0329 






0.0064 






0.0064 






0.0432 






(9 = 4) 


M 2 


0.0344 


1 


.2685 


0.0074 





.3890 


0.0073 





.4108 


0.0460 


1. 


.7936 




Mi 


0.0341 


1 


.3034 


0.0073 





.3824 


0.0072 





.4053 


0.0455 


1 


.8566 




B 


0.0432 


2 


.1693 


0.0105 





.5427 


0.0101 





.4901 


0.0546 


2 


.8808 


Gumbel 


Approx. 


0.0617 






0.0408 






0.0409 






0.0605 






(9 = 1.5) 


M 2 


0.0649 


2 


.2833 


0.0432 


2 


.1528 


0.0433 


1 


.9800 


0.0646 


3 


.0536 




Mi 


0.0639 


2 


.4663 


0.0426 


2 


.1943 


0.0428 


2 


.0367 


0.0638 


3 


.2590 




B 


0.0703 


2 


.5133 


0.0468 


1 


.8163 


0.0467 


1 


.7087 


0.0693 


2. 


.5655 


Gumbel 


Approx. 


0.0385 






0.0061 






0.0061 






0.0345 






(9 = 3) 


M 2 


0.0430 


1 


.9156 


0.0079 





.4884 


0.0074 





.3950 


0.0405 


2 


.5846 




Mi 


0.0426 


1 


.9784 


0.0078 





.4806 


0.0073 





.3858 


0.0400 


2. 


.6083 






0.0492 


2 


.2152 


0.0102 





.5375 


0.0099 





.4832 


0.0455 


2 


.1534 
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Table 3: Mean and MSE (xlO 4 ) Monte Carlo results. AR(1) and GARCH(1, 1) settings, 
sample size n = 200 and 1, 000 Monte Carlo replications. For each replication, we perform S = 2, 000 
tapered block multiplier (Mi) repetitions with Normal multiplier random variables, kernel function 
Ki, i = 1, 2, block length Im = 4, and block bootstrap (B) repetitions with block length Ib = 7. 



(ui,u 2 ) 




(1/3,1/3) 
Mean 


MSE 


(1/3,2/3) 
Mean 


MSE 


(2/3,1/3) 
Mean 


MSE 


(2/3,2/3) 
Mean 


MSE 


AR(1) setting with [3 = 0.25 






















Clayton 


Approx. 


0.0506 






n o^o 
u.uoou 






o n^n 

U.UOOU 






o n^n 
u.uoou 






(9 = 1) 


M 2 


0.0521 


1 


.5319 


U.UOOU 


1 
1 


.0527 


U.UOU1 


1 

1 . 


.1258 


o 0^4^ 
u.uo^o 


1 
1 


.4689 




Mi 


0.0517 


1 
1 


.6574 


0.0357 


1 


.1003 


0.0358 


1 


1662 


0.0541 


1 


.6083 




B 


0.0591 


n 
Z 


.0394 


0.0415 


1 


.4084 


0.0417 


1. 


.5256 


0.0624 


2 


.1027 


Clayton 


Approx. 


0.0272 






n 0048 
u.uu^to 






o nozis 

U.UU^to 






u.uoyo 






(9 = 4) 


M 2 


0.0285 





.8366 


n 00^9 

U.UUOZ 


u 


.2290 


U.UUOZ 


u. 


.2208 


041 3 
u.u^io 


1 


.3315 




Mi 


0.0283 





.8792 


n oo^i 

U.UUOl 


u 


.2273 


U.UUOZ 


u. 


.2200 


041 
U.U^AIU 


1 


.3858 




B 


0.0375 


1 


.9937 


0.0084 





.3726 


0.0085 


0. 


.3646 


0.0495 


2 


.1996 


Gumbel 


Approx. 


0.0518 






u.uo^ty 






u.uo^to 






0^07 
U.UOU I 






(9 = 1.5) 


M 2 


0.0549 


1 


.4150 


U.UO ( o 


1 
1 


.1973 


0377 
U.UO / / 


1 

1 . 


.2459 


0^0 

u.uoou 


o 

z 


.0067 




Mi 


0.0545 


1 


.5154 


n 0^70 

U.UO ( u 


1 
1 


.2339 


o n c t7' : t 

U.UO / 


1 
1 


.2871 


0^4^ 
U.UOOO 


o 
z 


.0766 




B 


0.0608 


o 
I 


.0500 


0.0415 


1 


.3234 


0.0418 


1 


,4093 


0.0608 


2 


.3118 


Gumbel 


Approx. 


0.0346 






n oo^& 

U.UUOo 






U.UUOo 






o n^o/i 
u.uou^ 






(9 = 3) 


M 2 


0.0386 


1 


.1686 


n nn7n 

U.UU ( u 


n 
u 


.3079 


U.UUOo 


n 

u. 


.2979 


u.uoou 


1 


.5135 




Mi 


0.0384 


1 


.2224 


n 0070 

U.UU ( u 


u 


.3052 


o nn^7 
u.uuo / 


u. 


.2892 


0^/17 
U.UO'i ( 


1 


.5157 




B 


0.0447 


1 

1 


.8434 


0.0095 





.3873 


0.0097 


0. 


.4303 


0.0409 


1 


.8797 


GARCH(1, 1) setting 
























Clayton 


Approx. 


0.0479 






U.UOOU 






u.uoou 






0^1 P. 
U.U01U 






(9 = 1) 


M 2 


0.0491 


1 


.1144 


n 0^47 


u 


.8485 


u.uo^o 


u. 


.8279 


0^90 

u.uozu 


I 

i 


.0958 




Mi 


0.0486 


1 


.3579 


O flQQO 

U.Uooy 


U 


.9021 


U.Uooo 


U 


.8100 


U.Uolo 


1 


.2013 




B 


0.0567 


2 


.0156 


0.0403 


1 


.1765 


0.0403 


1. 


,2556 


0.0600 


1 


.8542 


Clayton 


Approx. 


0.0252 






0.0055 






0.0056 






0.0403 






(9 = 4) 


M 2 


0.0259 





.5054 


0.0051 





.1979 


0.0053 


0. 


,2301 


0.0399 


1 


.0431 




Mi 


0.0258 





.6429 


0.0052 





.2199 


0.0051 


0. 


,2217 


0.0390 


1 


.0959 




B 


0.0345 


1 


.5252 


0.0081 





.2764 


0.0081 


0. 


.2921 


0.0484 


1 


.8359 


Gumbel 


Approx. 


0.0500 






0.0339 






0.0339 






0.0482 






(9 = 1.5) 


M 2 


0.0516 


1 


.0480 


0.0356 





.8486 


0.0354 


0. 


,8175 


0.0511 


1 


.2451 




Mi 


0.0516 


1 


.2235 


0.0351 





.9582 


0.0352 





.8848 


0.0503 


1 


.2774 




B 


0.0575 


1 


.5928 


0.0402 


1 


.2395 


0.0403 


1 


,2346 


0.0574 


1 


.9198 


Gumbel 


Approx. 


0.0341 






0.0074 






0.0074 






0.0291 






(9 = 3) 


M 2 


0.0362 





.9284 


0.0073 





.2941 


0.0072 


0. 


.2587 


0.0321 





.9133 




Mi 


0.0366 





.9782 


0.0073 





.2786 


0.0071 


0. 


,2888 


0.0320 





.9819 




B 


0.0435 


1 


.6811 


0.0103 





.3607 


0.0101 


0. 


,3390 


0.0390 


1 


.6878 
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respectively, and we chose the parameters in such a way that Kendall's r is either 
1/3 or 2/3, i.e., 6 6 {1,4} for the Clayton and 6 G {1.5,3} for the Gumbel-Hougaard 
copula. 

The AR(1) process we consider is the stationary solution of the AR(1) equation 

Xj = /3Xj_i + £j, 

where the innovations are supposed to have standard normal marginals linked by one 
of the aforementioned copulas and where the coefficient of the lagged variable is either 
j3 = 0.25 or (3 = 0.5. This stationary solution can be written as Xj = X^fcLo P ke j-k- 
We simulate an (approximate) sample of length n of this model as follows: for 
some reasonably large negative number M, e.g., M = —100, let JJj = (Uj t i,Uj t2 ), 
j = M, . . . , n be a sample of independent realizations of one of the aforementioned 
copulas. Set Bj = (*~ 1 (?7 fjl ) J *- 1 (t7 J - j2 )), with $ being the standard normal cdf, and 
recursively define Xjm = em and 

Xj = pX.j-i + Bj for all j = M + 1, . . . , n. (14) 

The last n observations form the sample Xi, . . . , X n . 

The GARCH(1,1) sample is simulated as following: define Bj = (sj i, £.7,2) as hi the 
AR(l)-example and recursively define om,% = \J^il{\ ~ a i ~ A) f° r * = 1, 2 and 

X hi = a j,i £ j,i' witn a j,i = ^i + Pirf-x,i + Q ' £ i-l,i ( 15 ) 

for i = 1,2 and j = M + 1, . . . ,n, where u x = 0.012, uj 2 = 0.037, A = 0.919, @ 2 = 
0.868, ai = 0.072 and a 2 = 0.115. Again, the last n observations form the sample 



Xi, . . . , X n . The considered coefficients are estimates derived in Jondeau et al. [271 ] to 
model volatility of S&P 500 and DAX daily (log-)returns in an empirical application 
which shows the practical relevance of this specific parameter choice. 

The theoretical covariance of Gc is easily calculated in the i.i.d. setting, but hardly deriv- 
able in the serially dependent case. Note that even a closed form expression for the copula 
of Xj is unaccessible. For that reason, we approximate the theoretical covariance by means 

of empirical covariances calculated from simulated samples of the process \fn{C n — Cat), 
where n = n(N) — > oo with n = o(N). We chose = 10 6 and n = 1000, and made 10 6 
replications on basis of which we calculated the empirical covariances. The results in Tables 
Q]-[3] (the 'True' vs. 'Approx.' lines) for the i.i.d. setting show that this approximation 
works sufficiently well, whence we can use it as a benchmark in the two serial dependent 
settings. 

Regarding the tapered block multiplier bootstrap we decided to use both kernel functions 
K\ and k 2 from Examples [3] and HI For the (moving) block bootstrap we choose the block 
length as Z#(100) = 5 and ^(200) = 7; this choice corresponds to /_e( n ) = |l-25n 1//3 J which 
satisfies the assumptions of the asymptotic theory. For a detailed discussion on the block 
length of the block bootstrap, we refer to Kiinsch (3o| as well as Buhlmann and Kiinsch 
0| . The tapered block multiplier technique is assessed based on a sequence (wj)j^z of 
normal random variables as introduced in Examples [3] and 2J The block length is set to 
l M (n) = [l.ln 1 / 4 ], hence Z M (100) = 3 and / A /(200) = 4, meaning that both methods yield 
2/jvf-dependent blocks. For all methods S = 2,000 bootstrap repetitions are performed and 
the target covariance is estimated by the sample covariance over the 2,000 repetitions. This 
procedure is repeated 1, 000 times and we report mean and mean squared error (MSE) for 
each method. 
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3. Testing for a constant copula 



The present section presents two nonparametric tests for a constant copula in the case of 
serially dependent processes. We begin with a test where a change point candidate is given, 
and proceed with a more general test without specifying a time point where a break in the 
copula structure occurs. Both tests are consistent against general alternatives and their 
finite sample performance is investigated by means of a simulation study. 



3.1. Specified change point candidate 

The specification of a change point candidate can for instance have an economic motivation: 
Patton 32fl investigates a change in parameters of the dependence structure between various 
exchange rates following the introduction of the euro on the 1st of January 1999. Focusing 
on stock returns, multivariate association between major S&P global sector indices before 
and after the bankruptcy of Lehman Brothers Inc. on 15th of September 2008 is assessed 



in Gaifier et al. 2l| and Ruppert [43|. Whereas these references investigate change points 
in functionals of the copula, the copula itself is in the focus of this study. This approach 
permits to analyze changes in the structure of association even if a functional thereof, such 
as a measure of multivariate association, is invariant. 

Suppose we observe a sample Xi, . . . , X n of a process (Xj)jgz. We derive a test for con- 
stancy of the copula in the case of a specified change point candidate indexed by [AnJ for 
A 6 (0, 1), i.e., a test for the hypothesis 



Hq : Uj ~ C\ for all j = 1, ... , n, 

C\ for all j = 1, . . . , [AnJ , 
C2 for all j = [\n\ + 1, . . . , n, 



where C\ and C2 are assumed to be different in at least one point u E [0, l] d (and hence, for 
continuity reasons, also in a neighbourhood of this point). The proposed test statistic will be 
based on a splitting of the sample into two subsamples: Xi, . . . , X^ ra j and X^ n j +1 , . . . , X n . 
A significant discrepancy between estimates of the copula in the two subsamples suggests 
to reject the null hypothesis. Assuming constant marginal distributions in each subsample, 
let 

^ 1 L An J ^ I n 

£l,...,LAnj(u)= p^T I^feup AAnj+l,...,nH = £ l {9j < u} 

1 J i=i L J i=LAnJ+l 

denote the corresponding empirical copulas of the two subsamples, where Uj.,... , Ur\ n J 

and V^ n j +1 , . . . , V n denote the pseudo-observations calculated from the first and second 
subsample, respectively. We consider the test statistic defined by 



T n (X) 



I UEM fr LX „ j( u)-e LA „ J+ ,..,„(u)} 



2 



du (16) 



which can be calculated explicitly; for details we refer to Remillard and Scaillet [40j|. These 
authors introduce a test for equality between two copulas which is applicable in the case 
of serial independence. Weak convergence of T n (A) under strong mixing follows from the 
following result. 
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Theorem 4. Consider observations Xi, . . . ,X n , drawn from a process (Xj)j g ^ satisfying 
the strong mixing condition o?x( r ) = 0{r~ a ) for some a > 1. Further assume a specified 
change point candidate indexed by [Xn\ for A E (0,1) such that Uj ~ C\, Xji ~ for 
all j = 1, . . . , [AnJ , i = l,...,d and Uj ~ C%, Xjj ~ i*2,i for all j = [Xn\ + 1, . . . , n, 
i = l,...,d. Suppose that C\ and C2 satisfy Condition {3)]. Under the null hypothesis 
C\ = C2, in the metric space £°°([0, l] d ), 

Y |Cl,...,LAnJ ~~ C^Xn\+l,...,n 

where Gc p (u) = Bc p (u) — Y2i=i D%C '(u)B(; (u^) forp = 1, 2 and where B^ andB<7 2 denote 
two independent tight centered Gaussian processes with covariances as specified in (|5|) . 



As a consequence of this Theorem, the continuous mapping Theorem yields 

T n (A) ^ T(A) = / { VT^AG C1 (u) - \/AG Ca (u) ) * du 
J[o,i] d L J 

for all A E (0, 1). Note that if there exists a subset X E [0, l] d such that 



y {vT^ACi(u) - \/AC 2 (u)} 2 du > 0, 



then T n (A) — > 00 in probability under H\. To estimate p- values of the test statistic, we 
use the tapered block multiplier technique developed in Section [2j For that purpose, let 
(^ni)j'eZ De a sequence of multipliers satisfying Al, A2, A3 and set 



where £n and £^ denote the arithmetic mean of the £j n in the corresponding samples. 

Let DjC n 1 and D^C n 2 denote estimators for the corresponding partial derivatives and set, 
for p = 1,2, 

d 

G£f n » = Bg nip (u) - £ A^(u)I^ n)P (u«). 



8=1 



In the i.i.d. multiplier case this statistic is, up to some negligible constants, the same as the 
one used in (4o| to test for equality between two copulas. 



Proposition 1. Consider observations Xj, . . . , Xu n j and X^ n j +1 , . . . , X„ drawn from a 
process (X,-),- e z- Assume that the process satisfies the strong mixing assumptions of Theo- 
rem^ Let [AnJ for A E (0, 1) denote a specific change point candidate such that Uj ~ C\, 
Xj_i ~ F\i for all j = l,...,|_An_|, i = l,...,d and Uj ~ C2, Xj^ ~ i 7 ^ for all 
j = [AnJ + 1, . . . , n, i = 1, . . . ,d. Suppose that C\ and C2 satisfy Condition |3|) and that 

D%C n \ and DiC n> 2 are estimators for the corresponding partial derivatives satisfying Cland 
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C2. Let £ n> i, . . . , £ n ,n denote samples of a tapered block multiplier process (£ n ,j)jeZ satis- 
fying Al, A2, A3 with block length l(n) — > oo, u>/iere Z(n) = C(n 1 / 2_e ) /or < e < 1/2. 
T/ien 

:= VT^XG^ - VXG^ A ^/T^AG Cl - v/AG C2 (17) 
in £°°([0, l] d ) 6oi/i under the null hypothesis as well as under the alternative. 



As a consequence, by the continuous mapping Theorem for the bootstrap, see 29j, 

ff(A)= / ft(u)) 2 rfuAT(A). 

J[0,l] d 5 

The integral involved in T* f (A) can be calculated explicitly [see [3£], Appendix B]. If we 

repeat the procedure S times to obtain a sample T^f^\ . . . ,T?f , then an approximate 
p-value for the test for Hq is provided by 

1 5 

5S 1 {l^W(A)>Tn(A)}- (18) 

Hence, p-values can be estimated by counting the number of cases in which the simulated 
test statistic based on the tapered block multiplier method exceeds the observed one. 



Finite sample properties. Size and power of the test in finite samples are assessed in a 
simulation study. We consider bivariate samples of size n = 100 or n = 200 generated as 
in Section [2,3[ i.e., marginal i.i.d., AR(1) and GARCH(1,1) processes are either linked by 
a Clayton or a Gumbel-Hougaard copula. The change point after observation [AnJ = n/2 
only affects the parameter within each family: the copula C\ is parameterized such that 
Kendall's t x = 0.2, the copula C 2 such that Kendall's r 2 = 0.2, . . . , 0.9. A set of S = 2, 000 
normal tapered block multiplier processes is simulated, where the kernel function is chosen 
as k>2 as suggested by the simulation results in Section 12.31 and where Zm( 100) = 3 and 
Zm (200) = 4 are chosen for the block length. 

The results of 1, 000 MC replications are shown in Tables [4] and [5] for n = 100 and n = 200, 
respectively. The test based on the tapered block multiplier technique leads to a rejection 
quota under the null hypothesis which is close to the chosen theoretical asymptotic size of 5% 
in all considered settings. Comparing the results for n = 100 and n = 200, we observe that 
the approximation of the asymptotic size based on the tapered block multiplier improves in 
precision with increased sample size. The tapered block multiplier-based test also performs 
well under the alternative hypothesis and its power increases with the difference t 2 — T\ 
between the considered values for Kendall's r. The power of the test under the alternative 
hypothesis is best in the case of no serial dependence as is shown in Table [5) If serial 
dependence is present in the sample then more observations are required to reach the power 
of the test in the case of serially independent observations. For comparison, we also show the 
results if the test assuming independent observations (i.e., the test based on the multiplier 
technique with block length Z = 1) is erroneously applied to the simulated dependent 
observations. The effects of different types of dependent observations differ largely in the 
finite sample simulations considered: GARCH(1, 1) processes do not show strong impact, 
whereas AR(1) processes lead to considerable distortions, in particular regarding the size of 
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Table 4: Size and power of the test for a constant copula with a specified change point 
candidate. Results are based on 1,000 Monte Carlo replications, n — 100, S = 2,000 tapered 
block multiplier repetitions, kernel function K2, and asymptotic significance level a = 5%. 
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the test. Results indicate that the test overrejects if temporal dependence is not taken into 
account; the observed size of the test in these cases can be more than twice the specified 
asymptotic size. For comparison, results for n = 200 and kernel function K\ are shown in 
Ruppert [43]. The obtained results indicate that the uniform kernel function K\ leads to a 
more conservative testing procedure since the rejection quota is slightly higher, both under 
the null hypothesis as well as under the alternative. Due to the fact that the size of the test 
is approximated more accurately based on the kernel function K2, its use is recommended. 

3.2. The general case: unspecified change point candidate 

The assumption of a change point candidate at specified location is relaxed in the following. 
Intuitively, testing with unspecified change point candidate(s) is less restrictive but a trade- 
off is to be made: the tests introduced in this section neither require conditions on the partial 
derivatives of the underlying copula(s) nor the specification of change point candidate(s), yet 
they are based on the assumption of strictly stationary univariate processes, i.e., Xji ~ Fi 
for all j € 7L and i = 1, . . . , d. The motivation for this test setting is that only for a subset 
of the change points documented in empirical studies, a priori hypothesis such as triggering 



economic events can be found [see, e.g., |l5j. Even if a triggering event exists, its start 



(and end) often are subject to uncertainty: Rodriguez [42| studies changes in dependence 
structures of stock returns during periods of turmoil considering data framing the East Asian 
crisis in 1997 as well as the Mexican devaluation in 1994, where no change point candidate 
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Table 5: Size and power of the test for a constant copula with a specified change point 
candidate. Results are based on 1,000 Monte Carlo replications, n — 200, S = 2,000 tapered 
block multiplier repetitions, kernel function K2, and asymptotic significance level a = 5%. 



T"2 


0.2 


3 

yj .O 


4 


Fi 


fi 


7 




8 




9 


i.i.d. setting 






















Clayton I = 1 


0.047 


0.172 


0.524 


0.908 


0.993 


1.000 


1 


.000 


1 


.000 


I = 4 


0.063 


n 1 a a 
0.164 


n ceo 


n fine 
U.905 


n non 
U.989 


i nnn 
1.U0U 


1 


.000 


1 

1 


nnn 
.000 


Gumbel I = 1 


0.043 


0.162 


0.525 


0.877 


0.991 


1.000 


1 


.000 


1 


.000 


I = 4 


0.055 


0.169 


0.535 


0.895 


0.996 


1.000 


1 


.000 


1 


.000 


GARCH(1,1) setting 




















Clayton I = 1 


0.040 


0.169 


0.503 


0.894 


0.994 


1.000 


1 


.000 


1 


.000 


I = 4 


0.056 


U.lOl) 


U.D41 


n nriQ 




l.UUU 


i 

i 


nnn 
.1)1)1) 


1 

1 


nnn 

.uuu 


Gumbel I = 1 


0.046 


0.154 


0.498 


0.873 


0.992 


1.000 


1 


.000 


1 


.000 


1 = 4 


0.057 


0.174 


0.496 


0.899 


0.994 


1.000 


1 


.000 


1 


.000 


AR(1) setting with j3 


= 0.25 


















Clayton I = 1 


0.052 


0.180 


0.521 


0.866 


0.989 


1.000 


1 


.000 


1 


.000 


I = 4 


0.057 


0.149 


0.497 


0.867 


0.988 


1.000 


1 


.000 


1 


.000 


Gumbel I = 1 


0.047 


0.180 


0.515 


0.872 


0.989 


1.000 


1 


.000 


1 


.000 


1 = 4 


0.050 


0.136 


0.490 


0.855 


0.992 


1.000 


1 


.000 


1 


.000 


AR(1) setting with /3 


= 0.5 


















Clayton I = 1 


0.107 


0.237 


0.523 


0.813 


0.975 


0.999 


1 


.000 


1 


.000 


I = 4 


0.058 


0.137 


0.396 


0.748 


0.957 


0.998 


1 


.000 


1 


.000 


Gumbel I = 1 


0.122 


0.227 


0.499 


0.825 


0.979 


1.000 


1 


.000 


1 


.000 


I = 4 


0.059 


0.123 


0.401 


0.770 


0.958 


0.998 


1 


.000 


1 


.000 



is given a priori. These objects of investigation are well-suited for nonparametric methods 
which offer the important advantage that their results do not depend on model assumptions. 
For a general introduction to change point problems of this type, we refer to the monographs 
by Csorgo and Horvath [i3| and, with particular emphasis on nonparametric methods, to 
Brodsky and Darkhovsky [2j. 

Let Xi, . . . , X n denote a sample of a process (Xj)jgz with strictly stationary univariate 
margins, i.e., Xji ~ Fi for all j G Z and i = l,...,d. We establish tests for the null 
hypothesis of a constant copula versus the alternative that there exist P unspecified change 
points Ai < . . . < Xp £ [0, 1], formally 

Hq : Uj ~ C\ for all j = 1, . . . , n, 

H\ : there exist = Ao<Ai<...<Ap< Ap+i = 1 such that Uj ~ C p 
for all j = [Ap-inJ + 1, . . . , [X p n\ and p = 1, . . . , P + 1, 

where, under the alternative hypothesis, C\, . . . , Cp+i are assumed to be pairwise different 
in at least one point u G [0, l] d . Unlike in the previous section, we estimate the pseudo- 
observations Ui, . . . , U n based on the whole sample Xj, . . . , X n . The following test statis- 
tics are based on a comparison of the empirical distribution functions of the subsamples 
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Ui,...,U^ n j and U^ n j +1 , . . . , U n : 

s (CU )- mt^mi f 1 vi 1 v i~ 

Observe the similarity between S n and the integrand in equation (fl~6|) of the previous section: 
except for the difference regarding the calculation of pseudo-observations, the functionals 
only differ by the factor W [C n \ ( n ~ [C n \ ) / n which assigns less weight to change point 
candidates close to the sample's boundaries. Define Z n : = {1/n, . . . , (n — l)/n}. We consider 
three alternative test statistics which pick the most extreme realization within the set Z n 
of change point candidates: 

rl = max/ S n ((,u) 2 dd n (u), (19) 
T% = max < jnax § n (C,u) — jxrin S n (C,u)>, (20) 

= l,...,n I 




T n 3 = max { max |S n (C,u)| } , (21) 
(ez n [u6{u f } i= i„ 

which are the maximally selected Cramer-von Mises ( CvM), Kuiper (K), and Kolmogorov- 
Smirnov (KS) statistic, respectively. We refer to Horvath and Shao [25j for an investigation 
of these statistics in a univariate context based on independent and identically distributed 
observations; T 3 is investigated in Inoue 26j for general multivariate distribution functions 
under strong mixing conditions as well as in Remillard (38| with an application to the copula 
of GARCH residuals. 



Theorem 5. Consider a sample Xi, . . . ,X n of a strictly stationary process (Xj)j g ^ sat- 
isfying the strong mixing condition ax(r) = 0(r -4-d { 1+e }) for some < e < 1/4. Then 
under the null hypothesis, in (£°°([0, l] rf+1 ), || • ||oo) , 



S n (C, u) 8 C (C, u) = B C (C, u) - CB C (1, u), 
where Mq(C, u) denotes a (centered) C -Kief er process, with covariance structure 

,<v} 



Cov(B c (Ci, u),B c (C2, v)) = mm(Ci, C2) £ Cov (l {Uo < u} , 1 {U 



for all Ci;C2 £ [0, 1] and u, v G [0, l] d . This in particular implies weak convergence of the 
test statistics 7r, 7r, and T 3 under Hq: 



-,1 w- 



^2 ™v 



sup / {B c (C,u)-CB c (l,u)} 2 dC(u), 

o<c<i J[o,i\ d 



T 

1 1 



3 w - 



sup 

o<C<i 



sup 

o<C<i 



sup {B c (C,u)-CB c (l,u)}- inf {B c (C,u) 
ue[o,i] d ue[o,i] d 



sup |{B c (C,u)-CB c (l,u)}| 

ue[o,i] d 



;i,u)} 
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Similar calculations as in [25[ reveal that — > oo for i = 1,2,3 under H\, Hence, a test 
which rejects Hq for unlikely large values of is consistent against general alternatives. 
Approximate critical values of the tests can be derived from the tapered block multiplier 
technique. 

Proposition 2. Consider a sample Xi, . . . , X n of a process (Xj)j- e z which satisfies Xji ~ 
Fi for all j Eli and i = 1, . . . , d. Further assume the process to fulfill the strong mixing as- 
sumptions of Theorem^ Let . . . , £ n , n denote a sample of a tapered block multiplier pro- 
cess (£,j, n )j£Z satisfying Al, A2, A3 with block length l(n) — > oo, where l(n) = 0(n 1//2_e ) 

for < e < 1/2 and define S n (C, u) = §^(C, u) - u), where 




(C,u)=n-^\ -^-1 1 



{U 3 <u}- 



Then, under the null hypothesis, S n — > Sc, while under the alternative § n = Op(l). 



An application of the continuous mapping theorem proves consistency of the tapered block 
multiplier-based tests. The p- values of the test statistics are estimated as shown in Equation 

OB. 



For simplicity, the change point location is assessed under the assumption that there is at 
most one change point. In this case, the alternative hypothesis can as well be formulated: 



H lb : 3X e [0, 1] such that Uj ~ 



C\ for all j = 1, ... , [Xn\ , 
C2 for all j = [An] + 1, . . . , n, 



where C\ and C2 are assumed to differ on a non-empty subset of [0, l] d . An estimator for 
the location of the change point A^, i = 1,2,3, is obtained by replacing max functions by 
argmax functions in Equations (I19p . (120j) . and (121j) . For ease of exposition, the superindex 
i is dropped in the following if no explicit reference to the functional is required. Given a 
(not necessarily correct) change-point estimator A n , the empirical copula of Xi, . . . , X^ ^ 

is an estimator of the unknown mixture distribution given by 

= l{x n <x}Ci(") + l^A}^ 1 [ACi(u) + (X - A) C 2 (u)] (22) 



[for an analogous estimator related to general distribution functions, see|l0|. The latter 
coincides with G\ if and only if the change point is estimated correctly. On the other 
hand, the empirical copula of X^ nj+l'"" ^ s an es thnator of the unknown mixture 
distribution given by 

C^ 2 (u) =l {Xn < A} (l - K)- 1 [(A - A n )d(u) + (1 - A) C 2 (u)] (23) 

+ 1 {A„>A} C '2(u), 

for all u G [0, l] d . The latter coincides with C2 if and only if the change point is estimated 
correctly. Consistency of A n follows from consistency of the empirical copula and the fact 
that the difference of the two mixture distributions given in Equations (|22p and (|23p is 
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maximal in the case A ra = A. Bai [lfl iteratively applies the setting considered above to test 
for multiple breaks (one at a time), indicating a direction of future research to estimate 
locations of multiple change points in the dependence structure. 

Finite sample properties. Size and power of the tests for a constant copula are shown in 
Tables [6] and [7] for n = 400 and n = 800, respectively. The observations are either serially 
independent or from a strictly stationary AR(1) process where the univariate innovations 
are either linked by a Clayton or a Gumbel-Hougaard copula. We consider the alternative 
hypothesis Hib of at most one unspecified change point. If present, then this change point is 
located after observation |_AnJ = n/2 and only affects the parameter within the investigated 
Clayton or Gumbel-Hougaard families: the copula C\ is parameterized such that Kendall's 
•n = 0.2, the copula C 2 such that Kendall's r 2 € {0.2,0.6,0.9}. We consider S = 1,000 
(n = 400) or S = 500 (n = 800) tapered block multiplier simulations based on normal 
multiplier random variables with block length Z^- (400) = 5, Im (800) = 6 and kernel function 

In the case of i.i.d. observations, we observe that the tapered block multiplier works simi- 
larly well as the standard multiplier (i.e., Im = 1): the asymptotic size of the test, chosen to 
be 5%, is well approximated and its power increases in the difference r 2 — T\. The estimated 
location of the change point, A n , is close to its theoretical value. Moreover, its standard 
deviation <r(A n ) as well as its mean squared error MSE(A n ) are decreasing in the difference 
T2 — T\. In the case of serially dependent observations sampled from AR(1) processes with 
(3 = 0.25, we find that the observed size of the test strongly deviates from its nominal size 
(chosen to be 5%) if serial dependence is neglected and the block length Im = 1 is used: its 
estimates are reaching up to 18.7%. The test based on the tapered block multiplier with 
block length £^(400) = 5 yields rejection quotas which approximate the asymptotic size 
well in all settings considered. These results are strengthened in Table [7] which shows results 
of MC simulations for sample size n = 800 and block length Im = 6. The power improves 
considerably with the increased amount of observations and the change point location is 
well captured. Standard deviation and mean squared error of the estimated location of the 
change point, A n , decrease in the difference t% — T\. 

Comparing the tests based on statistics T^, T%, and T^, we find that the test based on 
the Kuiper-type statistic performs best. The results indicate that the nominal size is well 
approximated in finite samples and that the test is most powerful in many settings. Like- 
wise, with regard to the estimated location of the change point, the Kuiper-type statistic 
performs best in mean and in mean squared error. 

The introduced tests for a constant copula offer some connecting factors for further re- 
search. For instance, Inoue (26| investigates nonparametric change point tests for the joint 
distribution of strongly mixing random vectors and finds that the observed size of the test 
heavily depends on the choice of the block length I in the resampling procedure. For differ- 
ent types of serially dependent observations, e.g., AR(1) processes with higher coefficient 
for the lagged variable or GARCH(1, 1) processes, it is of interest to investigate the optimal 
choice of the block length for the tapered block multiplier-based test with unspecified change 
point candidate. Moreover, test statistics based on different functionals offer potential for 
improvements. For instance, the Cramer-von Mises functional introduced by Remillard and 
Scaillet (4oT | led to strong results in the case of a specified change point candidate. Though 
challenging from a computational point of view, an application of this functional to the case 
of unspecified change point candidate(s) is of interest as the functional yields very powerful 
tests. 
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Table 6: Size and power of tests for a constant copula with unspecified change point can- 
didate. Results are based on 1, 000 Monte Carlo replications, n — 400, 5= 1, 000, kernel function 
K2, and a = 5%; additionally, the estimated change point location A„, <j(\ n ), and MSE(X n ) x 10 2 
are reported. 







T2 


size/power 
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0.9 


An 

0.6 




0.9 




M 

0.6 


0.9 


MSE 
0.6 


0.9 
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I = 1 


CvM 
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0.406 
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0.511 


o 


506 


o 
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0.061 


0.871 


0.372 






K 


0.040 


0.495 


o 
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0.496 


o 


496 


o 


074 


0.048 


0.556 
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KS 


0.062 


0.409 


o 
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0.507 


o 
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o 


079 


0.056 


0.627 


0.319 




I = 5 
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0.035 


0.342 


o 
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0.519 


o 
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o 


084 


0.058 


0.735 


0.349 
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0.031 


0.375 


o 
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0.495 


o 


496 


o 
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0.049 


0.489 


0.246 






KS 


0.036 


0.337 





881 


0.506 
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0.055 


0.645 


0.306 


Gumbel 


I = 1 


CvM 


0.047 


0.456 
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0.509 





507 





083 


0.055 


0.694 


0.304 






K 


0.056 


0.474 
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0.494 





495 
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0.050 


0.548 


0.251 






KS 


0.052 


0.437 





910 


0.502 





504 





080 


0.054 


0.644 


0.291 




I = 5 


CvM 


0.043 


0.387 





880 


0.506 





507 





086 


0.055 


0.739 


0.310 






K 


0.026 


0.343 





974 


0.487 





496 





072 


0.046 


0.518 


0.217 
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0.041 


0.349 





884 


0.497 
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0.642 


0.312 


AR(1) setting 
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= 0.25 
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1.390 


0.566 
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0.139 
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0.490 
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084 


0.052 


0.707 


0.276 






KS 


0.177 


0.499 





913 


0.504 
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0.071 


0.969 


0.503 




I = 5 


CvM 


0.046 


0.243 
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1.073 


0.497 
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0.496 
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0.040 
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KS 
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504 





089 


0.067 


0.793 
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Table 7: Size and power of tests for a constant copula with unspecified change point 
candidate. Results are based on 1, 000 Monte Carlo replications, n — 800, S — 500, kernel function 

n 2 , and a — 5%; additionally, the estimated change point location A„, ct(A„), and MSE(X n ) x 10 2 
are reported. 
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0.054 
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087 






KS 


0.046 
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068 


0.038 


0.463 
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0.585 
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243 
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0.037 
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0.032 


0.372 
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KS 


0.051 


0.509 
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0.502 
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0.044 
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199 




I = 6 


CvM 
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0.721 
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0.505 





504 





069 


0.037 


0.477 





138 






K 


0.034 


0.775 


1 


000 
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496 





059 


0.029 


0.355 





086 
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0.045 
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998 
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0.039 


0.428 
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AR(1) setting 
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= 0.25 
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0.184 


0.712 





995 


0.519 





507 





097 


0.056 


0.973 





321 






K 


0.122 


0.914 


1 


000 


0.495 





497 





063 


0.034 


0.411 





122 






KS 


0.169 


0.756 





999 


0.512 





506 





084 


0.050 


0.714 





248 




I = 6 


CvM 


0.065 


0.488 





939 


0.506 





508 





089 


0.054 


0.803 





293 






K 


0.057 


0.724 


1 


000 


0.493 





496 





060 


0.033 


0.359 





113 






KS 


0.062 


0.521 





965 


0.500 





503 





072 


0.045 


0.527 





206 


Gumbel 


I = 1 


CvM 


0.207 
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0.507 





509 
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0.052 


0.841 
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0.141 


0.879 
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000 
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497 





070 


0.036 


0.496 
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KS 


0.182 


0.780 
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0.500 





505 





089 


0.047 


0.805 
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I = 6 


CvM 


0.052 


0.533 
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0.514 





508 





083 


0.052 


0.707 





273 






K 


0.046 


0.670 


1 


000 


0.495 





496 





065 


0.032 


0.426 





106 






KS 


0.053 


0.520 





974 


0.508 





505 





076 


0.048 


0.584 





232 
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4. Conclusion 



Consistent tests for constancy of the copula with specified or unspecified change point 
candidate are introduced. We observe a trade-off in assumptions required for the testing: 
if a change point candidate is specified, then the test is consistent whether or not there is a 
simultaneous change point in marginal distribution function(s). If change point candidate(s) 
are unspecified, then the assumption of strictly stationary marginal distribution functions 
is required and allows to drop continuity assumptions on the partial derivatives of the 
underlying copula(s). Tests are shown to behave well in size and power when applied to 
various types of dependent observations. P- Values of the tests are estimated using a tapered 
block multiplier technique which is based on serially dependent multiplier random variables; 
the latter is shown to perform better than the block bootstrap in mean and mean squared 
error when estimating the asymptotic covariance structure of the empirical copula process 
in various settings. 
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Appendix A. Proofs 

Proof of TheoremUl Since Fi(Xj^) < itj if and only if Gi(Ujj) < Ui, where Gi denotes the 
empirical distribution function of Ux i, • • • , U n> i, we can assume without loss of generality, 
that Xj = Uj ~ C. Weak convergence of ~B>c,n under the strong mixing condition ax( r ) = 
0(r~ a ) for some a > 1 is established in Rio [4l| . Thus, Condition 2.1 in Biicher and 
Volgushev 0] is satisfied and an application of the functional delta method and of Theorem 
2.4 in 0| yields 

V^{C n (C^, C-i) - C} ^ G C , 
where C n j denotes the i-th marginal of C n , i = l,...,d. The assertion follows from 
||G c>n - V^{C n (C-{, . . . , C-J) - C}|U = 0(n-V2). . 

Proof of Theorem^ Without loss of generality we may assume Xj = Uj. Let C n ^ denote 
the empirical distribution function of the sample 

Uff 1+ i, . . . ,U Hi +i b ,~Uh 2 +i, ■ ■ ■ ,~U H2 +i B , >%+i> • • ■ ,~U Hk +l B - 

Again, \\G^ n - ^{C^C^, . . . , C^ J - C n }\\oo = O^- 1 / 2 ), and the result follows 
from Corollary 2.11 in [6j]. ■ 
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Proof of Theorem Without loss of generality we may assume Xj = JJj . The process 
= ^i'-t |f l fU ^, - C„(u,l =^t(|f-l) { 1(Uj£u) - C(u)} 



is defined on a product space (Qi x ^2, J~i ® ^2, Fi C5P2), where Xj is defined on (Oi, F\, 
and £ nj - on (^2,^2,^2) f° r all j = 1, . . . , n. It follows from the results in Biihlmann 



i 

Section 3.3] and [501 , Section 1.5] that there exists a set A\ 6 T\ with Pi(Ai) = 1 such that 
Bgjwi, •) B c in £°°([0, l} d ) for all oj x G Ai, 

where — ^ denotes weak convergence with respect to the multipliers. By Theorem 1.5.7 
and its addendum in [50| . B^f n (a;i,-) is asymptotically uniformly equicontinuous for each 

such u\. As sup u . g [ 0i i] \F~(m) — Ui\ = o(l), Pi-almost surely, we obtain 

M^ n (u)=M^ n (F-\u 1 ),...,F- 1 (u d ))^M c (Px-a.s.). 

Note that 



l;V -n¥ IL. =0(n~ 1 / 2 max |W&|) = op 2 (l), 



where the last estimation follows from P2(|£n,i| > s) = o(x~ 2 ) and P2(max™ =1 |£ n ,j| > 

eVn^n) < ^2(£n < 1/2) + nP 2 (|£n,i| > e^fn/2). Thus B^ n B c , Pi-almost surely, and 

since B^f n is a ball-measurable element of the space of cadlag functions [D ([0, l] rf ), || • ||oo) C 

(^°°([0, l] rf ), || • ||oo)) a detailed look at Theorem 1.7.2 and the proof of Theorem 2.9.6 in (5o| 
allows to translate this result into the Hoffmann- J0rgensen weak convergence, i.e., 

ic!„Y Bc m(r°([o,i] d ),|H|oo). 



;g B (u) = ig n (u) -^ac(u)b^(u« ' 

1=1 



Finally, set 



Proceeding as in the proof of Proposition 3.2 in [47] we obtain ||G^ n — Gp„||oo = °p(1)i 
which yields the assertion of the Theorem after an application of Lemma B.l in [5j. ■ 

Proof of Theorem^ Without loss of generality we may assume Xj = XJj. We begin by 
proving the joint weak convergence 

(B C ,n,l,Bc,n, 2 ) (Bc,i,B C)2 ) (A.l) 

in f°([0, l] d ) x £°°([0, l] d ), where 

LAnj 



C,n,l(u) = - 7 p=E{ 1 {U J <u}-C(u)} 



1 n 
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Asymptotic tightness of (ficnl>^>Cn2) follows from Lemma 1.4.3 in [50fl and the fact that 

<„,l + Op(l), 



1 



1 - a 



1 - S 



where both processes on the right-hand side are asymptotically tight by Theorem [TJ It 
therefore suffices to show that 

(Bc,n,l(ui), • • • ,B C ,n,l(uk),Bc,n, 2 (vi), . . . , Bc,n, 2 (vi)) 

^> (Bg,i(ui), . . . .BciCut), Bc^Cvi), . . . ,B c>2 ( Vl )) 

in lR fc+i , see Problem 1.5.3 in [501 ] . For the sake of a clear exposition we only consider the 
case k = I = 1, the general case follows along similar lines. We have to prove that, for each 



(ci,c 2 ) G 



ciB c ,n,i(u) + c 2 Bc,n,2(v) cil c ,i(u) + c 2 B Cj2 (v). 



Setting Zj(u) = l{u <u} — C( u )) the left-hand side of the previous expression can be 
written as ]Cj=i Cn,jj where 



Cnj' 



' Z j( U ) 1 {j<[Xn}} + / 2 , > r ^j( V ) 1 {i>LAnJ}- 



n 



^Jn- [Xn\ 



The asserted weak convergence follows from Theorem 2.1 in [341 ] . if we prove the corre- 
sponding conditions of that Theorem. We have 



E 



LAnj 



[An 



E B[Z ?; (v)Z,(v)] 



*J=1 



+ 



n — I An I 

J i,j=LAnJ+l 

|AnJ n 

- 2ClC2 V V E[Z ?; (u)Z,(v)]. 

>/LAnj(n-LAnj)fe. =L t3 +1 



The first two sums converge to 0^7(11,11) and c 2 7(v,v), respectively. Suppose that [AnJ > 
n — [Xn\ , the opposite case is treated analogously. Then, a tedious calculation shows that 
the third sum in the previous identity equals 



2cic 2 



^[Xn\(n-[Xn\) 



n— [AnJ |_AtxJ — 1 

kr(k) + ( n 

fc=l fc=n-|An| + l 



n-1 



|Ara|)r(£;) + ^ {n - k)r{k) 

k=\\n\ 



where r(|i — j\) = E[Zj(u)Zj(v)]. The first sum in the curly bracket multiplied with n" 1 
converges to by dominated convergence. If we multiply the other two sums with n _1 , 
then we obtain tails of absolutely converging series, which also converge to 0. To conclude, 



E 



( E ^ n >3' 

i=i 



7(u,u) 
7( v > v ) 



c = Var(ciB c ,i(u) + c 2 B c (v)). 



Since |Cnj| = 0{n~ 1 / 2 ) uniformly in j, some easy calculations show that condition (2.1) 
and (2.2) in 34] are satisfied. Hence, by Theorem 2.1 in that reference, the asserted weak 
convergence follows. 
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Define 

LAnj 



= vm^ 1 ^"' ~ C(U) L 



An application of the functional delta method to Mc, n ,p with the mapping from Theorem 
2.4 in [6] (and the usual estimation of the remainder term descending from the fact that 
F n j(Xij) < Uj and Xij < F~Hv,j) are not equivalent) yields 

(Gc,n,i,Gc,„,2) (G Cl ,G C2 ) in*°°([0,l] d ) x l°°([0,l] d ). 
Hence, Theorem 0] follows from the continuous mapping theorem observing that, under Hq, 



[Xn\ (n - [AnJ ) r - - •> ln-[Xn\~ [Xn 



n 



Ci,... ,LAnJ - Cn-LAnJ, 



Proof of Proposition^ Without loss of generality we may assume Xj = TJj. We consider 
the processes 



*C,n,l( u ) = 7=f I Tffi - 1 I 1 {U J -<u} 



1 

y/lXn 


LAnj 




1 

Vl Xn - 


LAnj 




1 




n 

E 

= LAnj+l 




LAnj . 


1 




n 

E 

= LAnJ+l 


\/n - 


LAnj j 



.1 



. sn 



{V^u}-^)} 



which are defined on the product space (Oi x r^-Fi <8>.F 2 ,Pi <8>P2)j where Xj is defined on 
(fil , J-"i , Pi ) and £ nj - on (^2, ^2^2) for all j = 1, . . . ,n. For the subsequent investigations, 

we may replace the arithmetic means £^ by 1 = E£ n j. As in the proof of Theorem [3j 
the results in Biihlmann 0, Section 3.3] guarantee the existence of a set A\ G J-"i with 
Pi (At) = 1 such that both 

B^,l(wi, •) B C)1 , B# n , 2 (wi, •) ^> B Ci2 for all u x € A x . 

(use an index shift for Bc,n,2)) where denotes weak convergence with respect to 

the multipliers. In the following we show joint weak convergence with respects to the 
multipliers, Pi-almost surely. Joint asymptotic tightness follows from asymptotic tightness 
of the components. As in the proof of Theorem H] it remains to show 

(Bg^Ul), • • ,ClW ) B^ 2 (v 1 ) ) . . . ,Bgn,2(vi))K •) 

^> (B 0>1 (ui), . . . ,Bo,i(uk),Bo^(vi), • • • ,B C)2 (vi)) 
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in which can be done by Theorem 2.1 in 34]. Again, for the sake of a clear exposition, 

we only consider the case k = I = 1. We have to prove that, for each c = (c\, c 2 ) £ M 2 , 



ciB$ n>1 (u) + c 2 B^ n>2 (v) ^> ciB c ,i(u) + c 2 B c , 2 (v) (P x -a.s.). 
The left-hand side of this display can be written as Xw=i Cn,j, where 



U,j — /, , 1 1 {U J <u} 1 {j<[AnJ} + 



An easy calculation shows that, uniformly in u, v and uj\, 

1 



L{U i <v}J-{i>LAnj}- 



|E 2 [Bg nil (u)B^„ )2 (v) 



< 



LAnJ n 

via»j(»-i^)§X. |Cw< ^ 



c-;(n)-i 



< const— 



i(n) 



0(1). 



Therefore, after some similar calculations as in the proof of Theorem[H <r 2 = E 2 [(^™ =1 ( n ,j) 2 ] 

Var(ciBc i i(u) + c 2 Bc i2 (v)), Pi-almost surely. Since |Cn,j| = 0(|^ n j — 1| • n -1 / 2 ) uniformly 
in u, v and UJ\, condition (2.1) in 34| is satisfied Pi-almost surely. Finally, we check the 
Lindeberg condition (2.2) in [|34|], i.e., for all e > 0, 



1 n 



n j=1 

for n — > 00. The expression on the left-hand side of this display can be estimated uniformly 
in u, v and oj\ by 

const • E 2 [(£n,i - l) 2 l{|e„ il -i[> CO r ls t v ^}] < const ■ -E 2 [(£ n) i - l) 4 ], 

which yields the assertion by the assumption on the finite centered moments of £ nj -. To 
conclude, 

(Bg n)1 ,Bg nt2 ) (Bc,i,B C)2 ) (Pi-a.s.) 
and analogously to the proof of Theorem [3] this implies 

(Bg n ,i,Ig n , 2 ) ^> (Bc,i,B c , 2 ) (Px-a.s.). 
Thus, by the continuous mapping Theorem 



«n,i(u (i) ; 



Vl - AGc,i(u) - v / AG C)2 (u) (Pi-a.s.). 



Again, the expression on the left-hand side is a ball-measurable element of D([0, l] d ) which 
allows to translate the latter weak convergence into the Hoffmann- J0rgensen conditional 
weak convergence. Finally, replace the true partial derivatives by their estimators and note 
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that the difference is uniformly op(l), which yields the assertion after an application of 
Lemma B.l in [5j. ■ 

Proof of Theorem [5l Without loss of generality we may assume Xj = Uj . Under the 
conditions on the mixing rate, it follows from Theorem 2 in Philipp and Pinzur [36] that 
the process 

Bc,„(C,u) := ^(l {x .< u} -C(u)} 

i=i 

weakly converges in £°°([0, l] a!+1 ) to a C-Kiefer process Bc(£, u). Hence, Condition 3.1 
in Q is satisfied and an application of Corollary 3.3 in that reference yields that 

S n (C, u) = — {\x n <F-^{u 1 ),...,x jd <F-^u d )} ~ CniC-Hm), c£ («*))} 

V 3=1 

weakly converges to §c(C,u) = Bc?(C,u) - C®c(l,u) in £°°{[0,l\ d+1 ). Under the null 
hypothesis we can write 

LCnJ 

s. 



(c^) = 4E{ 1 { u j <u } -^h} 



n 

3=1 

and the assertion follows from ||S n — 

Snlloo = P {n- 1 / 2 ). m 

Proof of Proposition [B Without loss of generality we may assume Xj = Uj . Suppose the 
null hypothesis holds and define 

LCnJ / ^ \ , LC«J 



We may replace by 1. Similar as in the proofs of Proposition [T] and Theorem[3]we begin 

by showing that B^f n — ^ Be, Pi- almost surely, where — ^ denotes weak convergence 
with respect to the multipliers. Convergence of the finite dimensional distributions follows 
along similar lines as in the previous proofs; the details are omitted for the sake of brevity. 
Tightness follows by tightness of Mq (1, •) and by analogous arguments as in the proof of 



Theorem 2.12.1 in 



5Qfl- 



Thus, as in the proof of Theorem [3l B r — > Be; and hence S„ — > S by the Lipschitz- 



"C,n ~ ^<-> 



continuous mapping Theorem for the bootstrap, see Proposition 10.7 in [29]. Under the 
alternative, again replacing £[^ n j by 1, we have 

, /LCnJ \ 2 j In 

-E 2 Yl " 1) 1{X,<U} < - £(1 " k / l nHk/ln) = 0(1) 

\j=l J 71 k=l 



by dominated convergence, which yields ||B^ n (^, u)||oo = Op 2 (l) uniformly in (,u and u±. 
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